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SECTION  1 
INTRODUCTION 

Structural  sandwich  composites  are  commonly  employed  in 
the  design  of  structural  components  requiring  high  stif fness-to- 
weight  ratio  characteristics.  Such  materials  are  particularly 
attractive  for  aerospace  applications  in  which  high  performance 
and  maximum  payload  are  typically  sought. 

While  sandwich  materials  offer  attractive  structural 
characteristics  for  many  applications,  effective  design  of 
sandwich  construction  is  generally  difficult  due  to  the 
complexity  of  response  and  the  possibility  of  failure  in  modes 
which  are  unique  to  sandwich  materials.  In  order  to  make 
effective  use  of  these  advanced  material  concepts,  improved 
analytical  techniques  are  needed  for  use  in  design  applications. 
The  scope  of  analysis  must  be  broad,  including  static  and 
dynamic  response  for  most  applications,  buckling  analysis 
for  sandwich  panels  and  shells,  and  thermal  stress  calculations 
for  heat-shielding  and  similar  components. 

The  development  presented  in  this  report  focuses  upon  such 

problems  as  static  response,  thermal  stress  determination  and 

natural  frequency  calculations,  all  of  which  may  be  repeated  a 

number  of  times  during  preliminary  design  and  are  thus  most 

useful  when  analyzed  in  an  interactive  mode.  Buckling,  large 

deflection  response  and  plastic  analysis  of  sandwich  often 

require  a  rather  detailed  and  time-consuming  analysis;  such 

nonlinear  problems  are  discussed  in  a  companion  report1.  Similar 

interactive  techniques  for  buckling  analysis  and  design 

2  3 

optimization  are  also  available  '  . 

1.1  GENERAL  APPROACH 

The  analytical  techniques  reported  here  are  oriented  toward 
layered  sandwich  panels,  including  stiffening  elements  and  edge 


T 


members,  and  are  limited  to  linear  elastic  responses.  Due  to 
the  practical  requirements  of  considering  many  different  types 
of  boundary  conditions  and  stiffener  arrangements,  and  possibly 
of  considering  non-rectangular  or  curved  panel  geometries,  the 
finite  element  method  is  chosen  as  the  analytical  basis 
throughout.  With  the  finite  element  approach,  only  the  task 
of  data  generation  need  be  specialized  to  particular  classes 
of  geometries  for  interactive  analysis. 

The  particular  finite  elements  employed  are  derived  from 
the  same  principles  as  the  nonlinear  elements  developed  in 
Reference  1.  A  typical  sandwich  panel  is  discretized  using 
thin  shell-type  elements  for  the  face  sheet  layers,  and  three- 
dimensional,  shear-flexible  elements  in  the  core  layer.  The 
shell  (face  sheet)  elements  are  formulated  using  a  penalty 
function  approach  which  eliminates  the  need  for  slope  compatibility 
between  elements  and  permits  the  joining  of  shell  (face  sheet) 
and  solid  (core)  elements  directly  in  a  straightforward  manner. 
Stiffening  members  derived  from  an  approach  similar  to  the  shell 
elements  can  be  used  to  model  either  face  sheet  stiffeners  or 
full-depth  spars  and  edge  closeouts.  In  each  case,  displacement 
compatibility  is  achieved  with  the  more  complex  elements  used  to 
represent  the  layers  of  the  sandwich  panel  itself.  Consistent 
static  and  thermal  loadings  are  provided  for  each  of  these 
classes  of  elements  in  static  analysis,  as  well  as  consistent 
mass  formulations  for  natural  frequency  calculation. 

1.2  SCOPE  OF  ANALYSIS 

In  the  present  development,  the  linear  (small  displacement, 
elastic)  behavior  of  sandwich  materials  is  considered  for  the 
cases  of  static  loading  and  free  vibration.  For  the  case  of 
static  analysis,  both  mechanical  and  thermally-induced  forces 
may  be  treated.  Provisions  are  made  for  performing  analyses 
of  models  of  moderate  size;  when  larger  analyses  are  necessary, 
the  data  generated  in  the  interactive  mode  may  be  used  in 
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conjunction  with  the  larger  finite  element  program  described 
in  Reference  1. 

Due  to  limitations  in  the  data  generation  phase  of  the 

computer  program,  the  class  of  sandwich  panels  considered  is 

flat  and  rectangular,  with  arbitrary  boundary  conditions  at 

each  edge.  However,  the  data  generation  scheme  is  organized 

so  that  geometries  for  non-rectangular  or  curved  panels  may  be 

generated  through  minor  modification  of  the  existing  code 

(see  Section  4).  Furthermore,  the  input  information  accepted 

by  the  analysis  program  is  in  the  formats  described  in 

Reference  1,  so  that  data  also  may  be  generated  manually  and 

4 

the  use  of  existing  interactive  plotting  software  is  possible. 

Each  layer  of  a  sandwich  panel  is  modeled  explicitly  in 
the  present  analysis,  so  that  localized  deformations  and  thermal 
stresses  due  to  non-uniform  heating  can  be  represented  with 
high  accuracy.  The  use  of  a  layered  model  also  facilitates 
the  representation  of  face  sheet  and  full-depth  stiffeners 
or  panel  edge  members.  These  stiffening  elements  may  be  located 
arbitrarily  within  the  sandwich  panel. 
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SECTION  2 

THEORETICAL  FORMULATION 


The  theoretical  basis  of  each  class  of  finite  element 
(face  sheet,  core,  stiffener)  used  in  the  present  development 
is  described  in  the  following  paragraphs.  Further  details  of 
the  element  formulations  may  be  found  in  References  2,  5  and  6. 
In  each  case  the  starting  point  is  the  element  potential  energy 
n  ;  for  static  problems,  the  true  state  of  equilibrium  is 
defined  simply  by  the  condition 


6np  -  °- 


(1) 


In  dynamic  (natural  frequency)  analysis  Hamilton's  principle 


t„ 
(  2 


(t  -  np)dv  =  o 


(2) 


provides  the  governing  equations,  in  which  T  is  the  kinetic 
energy. 


2.1  THIN  SHELL/FACE  SHEET  ELEMENTS 

The  development  of  shell-type  finite  elements  suitable 
for  the  representation  of  sandwich  face  sheets  is  complicated  by 
the  compatibility  requirements  for  elements  derived  from  classical 
thin  plate  or  shell  theories.  The  rotational  degrees  of  freedom 
dictated  by  the  condition  of  slope  continuity  preclude  the 
direct  enforcement  of  compatibility  with  standard  isoparametric 
solid  elements.  If  the  layers  of  the  sandwich  are  to  be  modeled 
explicitly,  it  is  necessary  either  to  derive  plate/shell  elements 
using  only  displacement  degrees  of  freedom  at  the  bond  lines  or  to 
construct  higher-order  solid  elements  to  represent  the  core 
layers.  In  the  present  work,  the  former  approach  is  adopted. 

The  thin  plate  or  shell  approximation  is  obtained  herein  by 
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appropriate  specialization  of  three-dimensional  continuum 
theory,  as  described  in  References  5  and  6. 

The  undeformed  geometry  of  a  typical  element  of  a  shell  is 
pictured  in  Fiqure  1.  Local  coordinates  x,y  are  imbedded  in 
the  base  plane  of  the  element,  and  the  height  of  the  middle 
surface  above  the  base  plane  is  denoted  by  Z(x,y).  Distances 
away  from  the  midsurface  are  measured  by  a  third  independent 
coordinate,  £.  The  position  of  a  generic  point  in  the  shell 
is  then 


r  =  xi  +  yj  +  zk  +  cn  (3) 

where,  if  the  shell  is  locally  shallow,  the  unit  normal  vector 
n  is 


n  =  -z,xi  -  z,y5  +  £•  (4) 

The  assumption  of  local  shallowness  also  implies  that  the 
coordinate  system  (x,y,c)  is  approximately  Cartesian. 

If  the  same  element  of  shell  material  is  subject  to  a 
displacement  field 


u  =  ui  +  vj  +  wk,  (5) 

the  displaced  position  of  an  arbitrary  point  is 

^  =  (x  -  £z,  +  u) i  +  (y  -  tz.  ,  +  v) j 

a  y 

+  (z  +  c  +  w)  k  .  (6) 

By  comparing  Equations  3  and  6,  the  three  dimensional  strains 
can  be  identified,  since 

2e^jdx^dx^  =  <?R  •  3r  -  c!r*  3r.  (7) 
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The  resulting  expressions  become,  after  linearization 


exx  “  u'x(1  ~  !’Z'xx^  ^Z'xyv'x  +  Z'xw'x 
eyy  *  v-y(1  '  ^-yy1  ‘  <z'xyu'y  +  z-yw’y 
eU  *  w's  -  Z'XU'C  -  z*yv'c 

2eyc  *  V'C(1  -  sz’yy>  +  w'y 

-  z.xu,y  ♦  z.y(w,?  -  v,y)  -  a,xyu,c 

2exc  ■  u-c  (1  •  CZ'xx>  +  w’x 

'  z'yv'x  +  Z'x(w't  '  u'x>  -  CZ'xyv'C 

2exy  *  u'y(1  -  «'xx>  +  v'x(1  ‘  'z'yy> 

+  z,  W  +  z,  w,  -  £Z,  (u,  +  V,  )  .  (8) 

x  y  y  x  xy  x  y 

The  behavioral  constraints  leading  to  thin  plate  or  shell 
response  are  expressed  directly  in  the  form 


ext  *  eyc  =  V  =  0  ’ 


Furthermore,  the  condition  of  vanishing  normal  stress  through 
the  thickness  is  imposed  exactly,  and  the  remaining  stress-strain 
relation  is  one  of  plane  stress. 


E 


Equation  10  is  written  symbolically  in  the  form 

a  =  De  (11) 

in  which  case  the  potential  energy  of  a  shell  is  simply 


=  |  |  eTDe  dV  -  j  uTt 
V  3V 


t  dA. 


(12) 


Here  3V  represents  the  shell  boundary,  and  t  the  components  of 
the  applied  mechanical  forces. 

To  complete  the  theoretical  formulation  of  the  element, 
the  constraints  summarized  in  Equation  9  must  be  introduced. 
This  is  accomplished  by  introducing  the  penalty  functional. 


F 


M 

Z 

i=l 


3 

Z 

j=l 


wij*Je(riJ) 


(13) 


In  Equation  13,  r^j  and  represent  collocation  points  and 
weights,  which  may  be  chosen  independently  for  each  of  the 
transverse  strains.  The  weights  play  the  role  of  penalty 
factors,  and  as  W^-h *°°  the  transverse  strain  constraints  are 
enforced  exactly  at  each  of  the  points  r^.  In  practice,  the 
magnitudes  of  the  VT  ^  are  chosen  to  provide  adequate  satisfaction 
of  the  constraints  and  good  numerical  conditioning  simultaneously 
for  the  practical  range  of  element  length-to-thickness  ratios*5. 
The  element  basis  is  simply  the  augmented  functional 


in  which 


n'  =  |  |  eTDe  dV  -  |  uTt  dA 


3V 

+  i  Z  y?Wy. 

2  i=1U~U 


(14) 
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rix> 


e 


yc 


) 


(15) 


and  W  is  the  matrix  of  weighting  factors  W 

In  thermal  stress  analysis,  the  constitutive  relation 
(Equation  11)  is  replaced  by 


a  -  D(e  -  aT) 


which  leads  to  the  energy 


IT  =  |  |  eTDe  dV  -  |  eTDaT  dV 


M 


-  |  uTt  dA  +  |  E  Y?WY, 


3v 


o  1 1  1  *1 

Z  i=l  ~1~  1 


(16) 


(17) 


The  second  term  of  Equation  17  leads  to  the  consistent  definition 
of  the  thermal  "loading"  induced  by  a  temperature  change  T(x,y,C) 
from  the  temperature  in  the  reference  state. 

For  dynamic  calculations,  the  kinetic  energy  is  formed 
directly  in  three-dimensional  form. 


T 


* 

p (5  -  r ) • (5  -  r )  dV. 
v 


From  Equations  3  and  6, 


(18) 


if  *2  • 2  *2 

T  =  i  p(u*  +  V  +  w  )  dV.  (19) 

v 

Since  the  potential  and  kinetic  energies  formulated  above 
involve  derivatives  no  higher  than  first  order,  the  finite  element 
approximation  must  satisfy  the  continuity  of  displacements 
(but  not  bending  slopes)  between  adjacent  elements.  In  the 
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present  study,  the  simplest  permissible  form  of  approximation, 
based  upon  linear  polynomial  shape  functions,  is  used.  The 
resulting  element  is  shown  in  Figure  2.  Despite  the  simplicity 
of  the  approximation,  the  element  exhibits  high  accuracy  in  both 
displacement  and  stress  computations,  and  is  particularly 
straightforward  to  use  in  comparison  with  many  other  shell 
finite  elements.  Linear  and  nonlinear  analyses  which  demonstrate 
the  accuracy  and  economy  of  the  present  shell  element  formulation 
can  be  found  in  Reference  6. 
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Similarly,  the  kinetic  energy  is 


•2  »2 
+  v  +  w  ; 


dV. 


(25) 


The  above  energy  expressions  are  discretized  using  a  linear 
polynomial  displacement  field  in  isoparametric  coordinates,  for 
compatibility  with  the  face  sheet  elements.  The  resulting 
finite  element  has  eight  nodes  (Figure  3)  and  24  degrees  of 
freedom. 

In  the  approximation  of  bending  response,  the  above 
eight-node  element  is  notoriously  poor  when  the  potential  energy 
(Equation  24)  is  evaluated  exactly**.  However,  the  stabilizing 
effect  of  surrounding  face  sheet  and  stiffener  elements  or 
closeouts  permits  the  element  energy  (i.e.,  the  stiffness 
matrix)  to  be  computed  using  a  reduced  integration  scheme  (one 
point  Gaussian  quadrature).  This  single-point  integration, 
though  inexact,  is  sufficient  for  convergence7  and  permits  a 
much-improved  simulation  of  the  panel  bending  response  over  the 
conventional  eight-node  element. 


2.3  FACE  SHEET  AND  FULL-DEPTH  STIFFENERS 

Discrete  stiffener  elements  for  use  with  the  present  face 
sheet  and  core  layer  elements  are  obtained  using  a  penalty 
function  approach  similar  to  that  employed  in  Paragraph  2.1,  in 
two  dimensions  rather  than  three.  The  resulting  element  is 
displacement-compatible  with  the  remaining  elements,  represents 
axial  and  primary  bending  effects  well,  and  may  be  arbitrarily 
oriented  in  three-dimensional  space. 

Both  types  of  stiffeners  are  shown  in  Figure  4.  In  each 
case,  the  stiffener  is  considered  to  be  relatively  small  in  width, 
so  that  all  stresses  are  approximately  uniform  in  the  width 
direction.  This  assumption  causes  inplane  bending  effects  to  be 
ignored,  which  is  consistent  with  the  resolution  provided  by  the 
face  sheet  discretization.  Normal  stresses  in  the  width  direction 
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are  also  neglected;  thus,  the  stiffener  is  subject  to  a  plane 
state  of  stress  in  the  remaining  dimensions,  although  arbitrary 
warping  and  twisting  of  the  element  is  permitted. 

The  strain  energy  of  a  typical  element  involves  only  the 
longitudinal  strains,  which  give  rise  to  both  axial  and  bending 
forces  in  the  deformed  state.  Remaining  components  of  strain 
(transverse  shear  and  normal  strain  in  the  depth  direction)  are 
suppressed  through  the  use  of  penalty  function  constraints 
(see  Paragraph  2.1).  Denoting  the  longitudinal  direction  by  x, 
the  potential  energy  is  simply** 


n 


p 


Eex*dv 


T 

Eci  e  dV 
xx 

V 


-  |  uTtdA 
8v 


(26) 


where  E  is  the  elastic  modulus,  and  the  kinetic  energy  is 


T  as  1 
2 


p(u2  +  v2  +  w2 ) dV . 


To  enforce  the  remaining  constraints 


(27) 


e 


zz 


introduce  the  penalty  functional 


(28) 


p  -  5  ”Az 


\  Vzz^z ) 


(29) 


in  which  W  ,  W  are  the  constraint  weights  and  r  ,  r  are  the 
corresponding  collocation  points.  Analogous  to  the  procedure 
of  Paragraph  2.1,  the  functional  F  is  added  to  the  potential 
energy  (Equation  26)  to  form  the  constrained  energy  on  which 
the  element  stiffenss  matrix  is  based; 
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"p 


-  I  \  Eexxav  -  1 


EaTe  dV 
xx 


+  5  V» (?x'  +  5  Vzz <Jz> 


-  J  uTtdV.  (30) 

3v 

The  above  potential  and  kinetic  energies  are  last  in 
discrete  form  through  a  linear  approximation  of  the  displacements 
u  within  individual  elements.  The  resulting  finite  elements, 
shown  in  Figure  5,  have  four  nodes  and  twelve  degrees  of  freedom. 
For  this  order  of  discretization,  the  element  strain  energy  may 
be  evaluated  exactly  using  a  four-point  (2  by  2  Gauss)  integration 
rule.  The  transverse  strain  constraints  (Equation  28)  are  each 
applied  at  a  single  point  corresponding  to  the  element  centroid 
location. 
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SECTION  3 

NUMERICAL  SOLUTION  TECHNIQUES 


The  methods  used  in  solving  the  governing  matrix  equations 
for  a  complete  finite  element  model  are  described  in  this  Section. 
Two  general  types  of  solutions  are  presented:  linear  simultaneous 
equations,  for  static  analysis  with  mechanical  and/or  thermal 
loading;  and  generalized  eigenvalue  systems,  for  use  in  natural 
vibration  analysis. 

3.1  LINEAR  STATIC  ANALYSIS 

In  static  analysis,  minimization  of  the  discrete  potential 
energy  for  an  assemblage  of  elements  leads  to  the  matrix 
equation 


KU  =  T  (31) 

in  which  K  is  the  stiffness  matrix,  U  are  the  (unknown)  nodal 
displacements,  and  T  is  the  load  vector  including  both  mechanical 
and  thermal  generalized  forces.  A  direct  solution  is  obtained  by 
a  Gauss-Doolittle  factorization  of  K,  followed  by  forward  and 
backward  substitution.  The  decomposition  of  K  is  of  the  form 

K  =  LDLT.  (32) 


Here  L  is  a  unit  lower  triangular  matrix. 


(nonzero) , 

i>j 

‘ij  -  - 

1, 

i=j 

0, 

i<  j 

(33) 

and  D  is  a  strictly  diagonal 

matrix.  Next 

let 

Z  = 

T 

DL  U 

(34) 
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and  the  relation 


LZ  =  T  (35) 

can  be  solved  for  Z  by  a  forward  substitution.  Prom  Equation  34, 

LTU  =  D-1Z  (36) 

which  can  then  be  solved  for  U  in  a  backward  substitution  step. 

Having  solved  for  the  nodal  displacements  U  over  the  entire 
finite  element  model,  a  complete  vector  of  displacements  can  be 
extracted  for  each  individual  element.  The  element  displacement 
vector,  along  with  the  shape  functions  for  an  element,  completely 
describe  the  state  of  displacement  (and,  therefore,  strain  and 
stress)  within  each  element  of  the  model. 

3.2  NATURAL  FREQUENCY  SOLUTION 

In  natural  vibration  analysis,  the  discretization  of  both 
the  potential  and  kinetic  energies  and  the  assumption  of 
harmonic  motions  yield  the  governing  system  of  equations 

KU  =  u)2MU  .  (37) 

Here  K  and  M  are  the  stiffness  and  mass  matrices,  respectively,  U 
is  the  nodal  displacement  vector,  and  oj  is  the  circular  frequency 
of  vibration. 

It  should  be  noted  in  Equation  37  that  the  mode  shape  U 
describes  only  the  relative  displacements  of  the  mesh  points 
(aU,  where  a  is  a  constant,  is  a  solution  as  well).  If  matrices 
K  and  M  are  of  order  N,  then  there  exist  N-solutions  of  Equation 

37:  uk,  U^ ;  i=l,2, . ,N.  Usually  only  the  lowest  several 

frequencies  and  mode  shapes  are  of  interest,  and  thus  a  vector 
iteration  procedure  is  most  appropriate  for  solving  the  system. 

An  approach  based  upon  simultaneous  vector  iteration  is  used  herein. 
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Factoring  the  stiffness  matrix  by  Choleski's  method  gives 

K  =  LLT.  (38) 

Equation  37  can  now  be  restated  as  the  standard  eigenvalue  problem 

AV  =  XV  (39) 

in  which 

-1  -T 
A  =  L  ML 

T 

V  =  L  U 

and 

X  =  -j  .  (42) 

U) 

The  eigenvalue  solution  is  performed  using  a  set  of  m  trial 
vectors,  where  m  is  greater  than  the  number  of  frequencies 
desired  but  is  generally  much  less  than  N.  These  trial  vectors, 
which  are  arranged  in  the  partial  modal  matrix 

*  -  [V.  V,  ...  .  V  ]  (43) 

gradually  converge  to  the  first  in  eigenvectors  (mode  shapes) 
as  the  iteration  proceeds. 

The  central  step  in  the  simultaneous  iteration  is  the 
function  of  the  m  by  m  interaction  matrix 

B  =  4>TA<!>.  (44) 

When  the  iteration  has  converged,  B  is  a  diagonal  matrix  since 
all  of  the  vectors  are  orthogonal  with  respect  to  A.  Before 


(40) 

(41) 
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the  solution  has  converged,  an  approximate  eigenvalue  analysis 
of  B  is  performed  at  each  iteration  to  generate  improved 
estimates  of  the  eigenvlaues  and  to  modify  the  iteration  matrix 
Details  of  the  implementation  of  this  procedure  may  be  found  in 
Reference  8. 
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SECTION  4 

COMPUTER  PROGRAM  DESCRIPTION 

The  finite  element  discretization  and  numerical  solution 
procedures  summarized  in  Sections  2  and  3  are  implemented  in  a 
computer  program  which  permits  the  interactive  solution  of  static 
and  natural  frequency  problems  for  sandwich  composites.  Auxiliary 
programs  are  used  in  conjunction  with  the  primary  analysis 
program  to  perform  the  tasks  of  data  input  and  generation,  and 
of  presentation  of  results. 

4.1  PROGRAM  FUNCTIONS  AND  ORGANIZATION 

The  system  of  computer  programs  which  are  used  for  the 
interactive  analysis  of  sandwich  panels  are  pictured  in  Figure  6. 

A  data  generation  program  is  used  to  prepare  data  interactively 
for  use  in  the  main  analysis  program.  The  sandwich  analysis 
input  data  (file  INPDAT)  conforms  to  the  input  data  formats  of 
Reference  1,  which  in  turn  are  patterned  after  those  of  the 

4 

MAGNA  general  purpose  finite  element  program  .  Thus,  the  geometry 
plotting  capabilities  described  in  Reference  9  may  be  used  to 
obtain  line  drawings  of  the  generated  mesh  for  documentation 
purposes  when  required.  Once  the  requested  analysis  has  been 
performed,  a  summary  of  results  is  stored  on  disc  (file  POSTFL) 

,  Q 

for  deformed  geometry  plotting  ,  listing  by  separate  output 
formatting  program,  or  listing  on  the  line  printer. 

4.2  INTERACTIVE  OPERATION 

The  interactive  programs  described  here  are  of  three  types: 
data  input  and  generation;  static  or  natural  frequency  analysis; 
and  tabular  output.  Operating  instructions  for  compatible  pre- 
and  post-processing  graphics  programs  are  given  in  Reference  9. 

Data  generation  is  performed  using  a  minimal  amount  of 
input  from  the  user  concerning  the  panel  geometry,  properties, 
supports,  loading,  temperature  distribution  and  mesh  divisions. 

In  particular,  the  following  data  are  requested: 
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1.  ANALYSIS  TYPE:  static  or  natural  frequency  solution. 


2.  NUMBER  OF  MODES:  number  of  natural  frequencies  and 
mode  shapes  to  be  solved  in  vibration  analysis. 

3.  TYPE  OF  LOADING:  mechanical  or  thermal  loads  (or 
both)  in  static  analysis. 

4.  PANEL  DIMENSIONS  A,B:  planform  dimensions  of  the 
sandwich  panel  in  the  x,y  directions.  At  present, 
the  data  generation  program  is  restricted  to  flat, 
rectangular  panels. 

5.  LAYER  THICKNESS  TL,TC,TU:  thicknesses  of  lower  face 
sheet,  core  layer,  and  upper  face  sheet. 

6.  FACE  SHEET  PROPERTIES  E , NU , ALFA , RHO :  Young's  modulus, 
Poisson's  ratio,  coefficient  of  thermal  expansion  and 
density  (requested  for  each  face  sheet  in  turn). 

7.  CORE  PROPERTIES  GX , GY , ALF AX , ALFAY , RHO :  shear  moduli, 
coefficients  of  thermal  expansion  and  density  for 
the  core  layer. 

8.  NUMBER  OF  STIFFENERS  NSX,NSY:  number  of  stiffening 
members  in  the  x  and  y  directions,  respectively.  Both 
face  sheet  and  full  depth  stiffeners  are  included. 

9.  STIFFENER  LOCATIONS  X(OR  Y) ,  I TYPE:  for  each  stiffening 
member,  the  constant  value  of  x  or  y  at  which  the 
member  is  located  is  requested?  valid  x-values  are  between 
0  and  A,  and  valid  y-values  are  between  O  and  B.  For 
each  member,  the  type  is  also  requested  (lower  face, 
upper  face,  full  depth). 

10.  STIFFENER  PROPERTIES  E,  ALFA,  T,  RHO:  for  each 
stiffener,  these  are  the  values  of  the  Young's  modulus, 
the  coefficient  of  thermal  expansion,  the  stiffener 
width,  and  the  mass  density. 

11.  SUPPORT  CONDITIONS  ISO, ISA, JSO, JSB:  for  each  of  the 
edges  of  the  panel  (x  =  o,a;  y  =  o,b),  a  boundary 
condition  must  be  specified.  Support  conditions  which 
can  be  chosen  independently  at  each  panel  edge  include: 
free,  simple  support,  clamped,  symmetric  and 
antisymmetric . 

12.  PRESSURE  LOADINGS  PL,PU:  in  linear  static  analysis, 
pressures  may  be  specified  on  the  lower  or  upper  face 
sheet  surfaces,  or  both.  Note  that  all  pressure  values 
are  positive  upward. 
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13.  POINT  LOADINGS  XP,YP, FORCE:  location  and  magnitude 
of  point  loadings  for  static  analysis. 

14.  NUMBER  OF  ELEMENTS  PER  EDGE  NEX , NEY :  this  data 
defines  the  degree  of  mesh  refinement  in  the  model 
by  dividing  the  sandwich  panel  into  a  rectangular 
grid  of  NEX  by  NEY  finite  elements.  Usually  a 
moderate  number  of  elements  (4-6  per  edge)  provides 
a  good  balance  of  accuracy  and  solution  cost. 

15.  GRID  POINT  LOCATIONS  XG  (OR  YG) :  along  each  side  of 
the  panel,  the  locations  of  grid  lines  may  be  specified 
to  generate  nonuniform  (graded)  meshes.  NEX-^1  locations 
must  be  specified  along  the  x-axis,  and  NEY+1  locations 
on  the  y-axis.  Note  that  each  of  the  stiffener 
locations  specified  in  Item  9  must  correspond  to  a  grid 
line  location. 

In  addition  to  geometric  parameters  and  other  data  defined 
in  the  interactive  mode,  temperature  data  must  be  supplied  if 
thermal  stress  analysis  is  to  be  performed.  By  default,  the 
data  generator  will  create  temperature  data  for  a  uniformly 
heated  panel  (T  =  100°,  in  consistent  units).  To  specify  a 
different  temperature  distribution,  the  user  must  supply  a 
separate  subroutine  to  generate  the  required  data  T  =  T(x,y,z). 

This  user-written  subroutine  has  the  form  shown  in  the  following 
example,  for  which  the  temperature  distribution  is 


T (x,y,z ) 


1 . 2Ze 


-0.028x 


(45) 


The  appropriate  user  subroutine  (which  always  has  the  name 
UTEMPS)  is 


SUBROUTINE  UTEMPS  (X,Y,Z,T) 

T  =  1. 2*Z*EXP (-0. 028*X) 

RETURN 
END  . 

User  subroutine  UTEMPS  is  called  once  for  every  generated  point 
in  the  finite  element  model. 

The  analysis  is  next  performed  by  simply  loading  and 
executing  the  analysis  module.  While  the  final  results  of  an 
analysis  are  written  to  the  file  POSTFL  for  later  processing, 
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certain  information,  including  diagnostic  and  warning  messages, 
is  output  directly  at  the  user's  terminal.  In  the  event  that  a 
problem  becomes  too  large  for  the  interactive  program,  the 
analysis  can  be  performed  in  the  batch  mode  using  the  program 
of  Reference  1,  since  the  input  data  are  fully  compatible. 

Following  the  analysis  step,  tabular  results  may  be 
printed  by  simply  executing  the  output  module.  Analysis  results 
provided  in  table  form  include  deflection  data  on  selected  grid 
lines,  and  strain  and  stress  data  at  individual  element  centroids. 
Alternatively,  the  file  POSTFL  may  be  used  as  input  to  graphics 
programs  which  produce  deformed  geometry  plots  for  each 
visualization  of  the  problem  results. 
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SECTION  5 
SAMPLE  ANALYSES 


Typical  analyses  performed  using  the  solution  techniques 
described  in  Sections  2  and  3  are  presented  in  this  section. 

Static  problems  involving  mechanical  and  thermal  loading,  as 
well  as  natural  frequency  solutions,  are  discussed.  The  operation 
of  the  analysis  program  is  demonstrated  in  greater  detail  in  the 
sample  terminal  session  included  in  the  Appendix. 


5.1 


STATIC  ANALYSIS  WITH  MECHANICAL  AND  THERMAL  LOADS 


To  demonstrate  the  static  analysis  option,  a  rectangular 
panel  previously  studied  by  Monforton^®  and  Kan  and  Huang^  is 
considered.  The  flat  sandwich  plate  has  the  following  geometry 
and  properties: 


Planform: 


Thicknesses: 


Face  Material:  E. 


a=b  =  50  inches 

tf  =  0.015  inches 

t  «  1.000  inches 
c 

,6 


Core  Material: 


10.5  x  10°  lb/in 

vf  =  0.30 

a  £  =  12.4  x 

G_  =  50000. 
c 

a  =  12.4  x  10  in/in°F 


10  ®  in/in°F 


lb/  in 

.-6 


Using  these  parameters,  a  static  solution  is  performed  for 
the  case  of  a  fully  clamped  boundary  and  uniform  lateral  pressure 
(Figure  7) .  Due  to  symmetry  of  the  problem  only  one  quarter 
of  the  panel  is  modeled.  The  mesh  consists  of  a  5  by  5  grid 
of  elements  in  each  layer,  for  a  total  of  260  unconstrained 
degrees  of  freedom.  For  a  unit  lateral  pressure,  the  computed 
central  displacement  is  cjc  =  0.113  inches.  This  value  is  in  good 
agreement  with  previously  published  results;  for  example, 
linearization  of  the  results  of  Reference  11  gives  ojc  =  0.09497 
inches . 


27 


A  static  analysis  has  also  been  performed  for  the  above 
panel  considering  thermal  loading  only.  In  this  case,  the 
panel  is  assumed  to  be  simply  supported.  The  temperature 
distribution  is  such  that  the  two  face  sheets  are  at  opposite 
uniform  temperatures,  and  the  core  experiences  a  linear  thermal 
gradient.  The  following  user  subroutine  is  used  to  define  the 
thermal  field: 

SUBROUTINE  UTEMPS  (X,Y,Z,T) 

T  =  200. *Z 

IF  (ABS ( Z) • GT . 0 . 5) T  =  100 . *T/ABS (T) 

RETURN 
END  . 

This  antisymmetric  temperature  field  produces  a  pure  moment  type 
of  loading,  and  thus  a  bending  response.  The  resulting  deformed 
shape  of  the  panel  is  shown  in  Figure  8.  Typical  Stress  results 
(in  this  case  the  core  shear  stress  Txz)  are  depicted  in  graphical 
form  in  Figure  9. 

5.2  NATURAL  FREQUENCY  ANALYSIS 

A  simply  supported  sandwich  panel  is  considered  to  determine 
its  free  vibration  response.  The  following  geometric  and 
material  parameters  are  used: 

Planform:  a  =  62.25  inches 

b  =  43.50  inches 

Thicknesses:  t^  =  0.072  inches 

t  =  1.856  inches 

c 

Face  Sheets:  Ef  =  10.  x  106  lb/in2 
v  f  =  0.33 

=  0.000244  lb-sec2/in4 

Core:  Gc  =  30000.  lb/in2 

Pc  =  0.0000123  lb-sec2/in4  . 

The  finite  element  model  of  one  quarter  of  the  panel  is  shown  in 
Figure  10.  Thirty-six  elements  are  used  in  each  layer  of  the 
model. 
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Figure  10.  Sandwich  Panel  for  Natural  Frequency 
Analysis. 


The  solution  for  the  fundamental  frequency  of  the  panel 

using  the  model  shown  yields  w.  =  174  Hz.  This  value  compares 

1  12 

well  with  the  classical  analytical  solution  of  =  170  Hz. 

A  second  solution  has  also  been  performed  for  a  similar 
panel  of  larger  aspect  ratio.  Panel  data  for  this  example 
is  identical  to  the  problem  above,  with  the  following  exceptions 

b  =  34,00  inches 
t^  =  0.063  inches 

t  =  1.874  inches. 

In  this  case,  the  computed  frequency  of  240  Hz  is  also  in  good 

.  12 
agreement  with  a  closed-form  solution  of  =  250  Hz. 


SECTION  6 

SUMMARY  AND  CONCLUSIONS 


A  finite  element  approach  for  the  static  stress  analysis 
of  sandwich  panels  under  mechanical  and  thermal  loadings, 
and  for  the  natural  frequency  analysis  of  sandwich,  has  been 
described.  Layers  of  a  sandwich  panel  are  modeled  explicitly, 
so  that  unequal  face  sheets,  stiffening  members,  and  other 
design  details  are  represented  with  relatively  high  accuracy. 

Due  to  the  finite  element  basis  of  the  present  methodology, 
rather  arbitrary  geometries  and  variations  in  mechanical  pro¬ 
perties  are  easily  considered. 

The  methods  described  are  suitable  for  interactive  use, 
to  facilitate  parametric  and  design  studies.  Data  generation 
utilities,  which  permit  the  definition  of  simple  panel  geometries 
from  a  minimal  amount  of  input,  are  provided  to  support  this 
mode  of  usage.  Generated  input  is  compatible  with  that  of  a 
comparison  program1  with  larger  capability.  Results  are  out¬ 
put  in  a  form  suitable  for  plotting  using  existing  computer 
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graphics  programs. 

The  developments  described  herein  can  be  improved  to 
provide  more  detailed  and  comprehensive  analytical  capabilities 
(e.g. ,  curved  panels,  arbitrary  stiffener  shapes,  linear 
buckling  analysis,  etc.).  However,  in  their  present  form, 
the  methods  presented  represent  a  valuable  capability  for 
practical  design  analysis  of  sandwich  composite  structures. 
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APPENDIX 

SAMPLE  INTERACTIVE  TERMINAL  SESSION 

A  typical  interactive  analysis  session  using  the  sandwich 
analysis  program  described  herein  is  reproduced  in  this  Appendix. 
This  sample  session  demonstrates  the  procedures  for  accessing 
and  executing  the  data  generation,  analysis  and  plotting 
programs,  and  illustrates  the  types  of  output  which  can  be 
obtained  from  the  analysis. 

The  analysis  procedure  is  initiated  by  accessing  a  stored 
control  language  procedure  file  (denoted  by  "PROCFIL"  in  the 
example) ,  and  issuing  a  command  to  start  the  procedure: 

BEGIN, SWICH , PROFIL . 

This  command  automatically  invokes  the  data  generation,  analysis 
and  tabular  output  modules.  By  appending  the  keyword  "PLOT" 
to  the  BEGIN  command,  a  geometry  plotting  utility  is  also 
executed  before  and  after  the  analysis  module.  The  "PLOT" 
option  is  used  in  the  sample  session  to  clarify  the  problem 
geometry  and  results. 
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